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My life flows between confines 

But knowledge has no confines. 

If we use the constrained 

To follow after the unconstrained, 

There is a danger that the flow will cease; 

And when it ceases, 

To exercise knowledge is purest danger. 
Chuang-tzu 

1. INTRODUCTION 

The isocurvature baryon model is one of only a handful of explanations for structure formation in the 
universe (Suginohara & Suto 1992; Cen, Ostriker & Peebles 1993). This model was proposed by Peebles 
(1987) to satisfy dynamical measurements of a low fio ~ 0.1 — 0.2 (e.g. Peebles 1986, Dekel et al. 1993) while 
simultaneously avoiding the need for hypothetical non-baryonic dark matter or a non-zero cosmological 
constant. Recent observations of microlensing from compact halo objects may also support the idea of 
baryonic dark matter (Alcock et al. 1993). With the COBE DMR detection of large angle anisotropics in the 
Cosmic Microwave Background (CMB) and the COBE FIRAS upper limit on Compton-t/ spectral distortions, 
however, the isocurvature baryon model with the simplest and most plausible thermal histories is ruled out. 
Normalized to the COBE DMR detection at 10° of (AT/T) rms = 1.12 ± 0.10 x 10" 5 (Bennett et al. 1994), 
the isocurvature baryon model with a standard recombination history vastly overproduces temperature 
fluctuations just under the degree scale. In this model, however, structure can form immediately after 
recombination, and an early reionization is plausible (Peebles 1987). Nevertheless, although rescattering 
damps the small scale temperature fluctuations to allowable levels, it also gives rise to unacceptably high 
spectral distortions in a fully ionized model (Gnedin & Ostriker 1992, Tegmark & Silk 1994). Therefore, 
only models with a non-standard thermal history, e.g., a partially or recently ionized universe, can survive 
these opposing constraints. 

Moreover, it is by no means obvious that one can construct a model consistent with CMB observations 
at the degree scale even allowing for an ad hoc ionization history. Recently, Chiba, Sugiyama, & Suto (1994) 
have shown that the fully ionized isocurvature baryon model is inconsistent with the SP91 13 point scan 
detection at ~ 2° (Schuster et al. 1993). Indeed we show here that if low fluctuations at this scale are 
confirmed, all isocurvature baryon models are excluded regardless of ionization history or even the possible 
formation of compact objects such as black holes. However, several recent CMB experiments, e.g. Python 
(PYTH, Dragovan et al. 1993), MAX GUM (Gundersen et al. 1993), MSAM (Cheng et al. 1993), have all 
reported higher detections than SP91 on similar angular scales. Moreover, the second year of COBE DMR 
data and the high detection by Tenerife (TENE) at the 5° scale (Hancock et al. 1993) indicate a steeper than 
Harison-Zel'dovich spectral slope at large scales. The standard isocurvature baryon model, which has low 
fio, can naturally account for this effect (Sugiyama & Silk 1994). Furthermore, since large scale fluctuations 
are nearly independent of ionization history, all models considered here can also account for this feature. 

Since the data at scales smaller than the COBE DMR detection are marginally consistent at best, 
it is no surprise that no isocurvature baryon model can satisfy all the constraints from the data. Until 
this situation is resolved, we cannot judge the viability of this model on the basis of intermediate scale 
anisotropics. We do however present the full CMB predictions of these models which can be employed once 
the measurements improve. In light of this current uncertainty, we also place constraints on the thermal 
history of this model from first and second order temperature fluctuations at arcminute scales, spectral 
distortions, and galaxy clustering. We provide the first consistent treatment of the matter power spectrum 
in compact object dominated universes and find significant differences from other works (e.g. Cen, Ostriker, 
& Peebles 1993), which may have undesirable consequences for large scale structure formation. The favored 
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values for the slope of the power spectrum n > —0.5 also overproduce galaxy clustering in the absence of 
antibias, if consistency with spectral distortions is also required. 

The plan of this paper is as follows. In §2, we explore the dependence of matter and temperature 
fluctuations, as well as spectral distortions, on thermal history including compact object formation. This 
allows us to place constraints and make predictions for the various models in §3. In §4, we discuss the 
implications of these results and evaluate the present status of the isocurvature baryon scenario with respect 
to CMB anisotropics. 

2. THERMAL HISTORY EFFECTS 

First, we briefly summarize the general properties of the isocurvature baryon model. This model consists 
of only baryons, photons and three massless neutrinos in a low density open universe. Initially, there are no 
fluctuations in the total density, but rather perturbations in the entropy per baryon. We thus assume initial 
conditions such that this entropy fluctuation follows a pure power law in k, S 2 (k) oc k n where k 2 = k 2 + K 
and K = —Hq(1 — flo) (Wilson 1983). In such models, the absence of Silk damping on small scales allows 
for the collapse of objects immediately following standard recombination at z ps 1000. This first generation 
of objects could reionize the universe to a significant degree (Peebles 1987). It is also possible that a 
large fraction of the mass ends up in these compact objects and thereupon behaves effectively as cold dark 
matter (Gnedin & Ostriker 1992). Gnedin & Ostriker also propose that compact object formation may 
alter nucleosynthesis and resolve the most serious difficulty of this model: fio = &b ^ 0.015/} -2 , the value 
required by standard nucleosynthesis (Smith et al. 1993). A second phase of structure formation, perhaps 
associated with galaxy formation, could further ionize the universe. Therefore the following parameters 
adequately describe the range of possibilities for the model: z c , the epoch of compact object formation; 
x c , the ionization fraction after z c ; £Iigm, the amount of matter left in the intergalactic medium; z, the 
secondary ionization redshift; and x,, the ionization fraction after secondary ionization. We will now describe 
the effects that these parameters have on the evolution of perturbations in the matter and radiation. For a 
more detailed treatment of these effects, see Hu & Sugiyama (1994). 

2.1 The Matter Power Spectrum 

Let us first consider the simplest case where the effect of compact objects is negligible and reduce our 
parameters to x, and z,. Gravitational instability and microphysical processes will transform the initial 
perturbation into a power spectrum of the form P(k) = (T(k)S(k)) 2 . We solve the first order Boltzmann 
equations for the coupled photon-baryon system as well as the neutrinos, to the present, in order to obtain 
the transfer function T(k) (Sugiyama & Gouda 1992, Hu & Sugiyama 1994). 

Even with isocurvature initial conditions, a curvature perturbation will be stimulated in order to keep 
the entropy constant (Kodama & Sasaki 1986). These perturbations will grow as the adiabatic mode and 
become important as the perturbation enters the horizon in the matter dominated epoch. However, when 
the adiabatic component enters the Jeans length, pressure gradients will force it to oscillate and damp away. 
Since entropy is conserved in the tight coupling limit, on small scales the baryons are left with the initial 
entropy perturbation. On large scales, the initial entropy perturbation is shifted from the matter to the 
radiation as the universe becomes matter dominated, in order to avoid a large curvature perturbation on 
superhorizon scales. Thus the matter is left with a characteristic k 2 (l — 3K/k 2 ) tail to the transfer function 
due to the remaining fluctuations from the adiabatic growth. 

The peak of the transfer function therefore corresponds to the maximal Jeans scale (see Fig. 1). For 
models in which standard recombination is followed by a significant transparent period, the peak will be close 
to the Jeans scale at recombination, k ( j SR) k, 0.13[1 + 0.24(fi /j 2 ) _1 + 0.007(fi o fc 2 ) _2 ] 1/2n ofc 2 Mpc _1 , and the 
oscillations at small scales will not have had a chance to damp away. For a universe that never recombines 
(or z, ?s 1000), the maximum Jeans length is larger, k^j NR ^ ps 0.13f2o^ 2 Mpc _1 , since standard recombination 
is close to matter-radiation equality in low fio models. Between standard recombination and z,, fluctuations 
grow in linear theory, and thus will be larger for models with low z,. On the other hand, from z, to the end 
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of the drag epoch Zd rag ps 160(f2o^ 2 ) 1 ^ 5 *j- 2 ^ 5 , Compton drag will again prevent fluctuations from growing 
inside the Jeans length. Since Compton drag is less effective for low x,, these models have larger fluctuations 
at the present and a less prominent peak. If claims that the observational power spectrum is quite smooth 
are confirmed (Peacock & Dodds 1994), then these models may be the only viable ones. 
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FIG. 1. Matter transfer function in partially ionized models. Scenarios with low z, have their peak on 
smaller scales due to the smaller Jeans length at standard recombination. Those with high x, suffer 
suppression in growth due to high Compton drag. Compact object dominated scenarios will resemble 
the standard recombination model of z, =0. 
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FIG. 2. Present day amplitude of fluctuations on the 8/i _1 Mpc scale <r%(n) for a normalization to the COBE 
DMR detection at 10°. Fluctuations in the no recombination (NR, thick) scenario are suppressed with 
respect to the standard recombination counterparts (SR, thin) due to Compton drag. The value of fio 
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decreases from top to bottom for h = 1.0 (solid line) and h = 0.5 (dashed line). All partially ionized 
scenarios have fluctuations in between these limits. Models which are compact object dominated have 
nearly the (SR) values. Analytic fits to these numerical results are found in Tab. 1. 

The present amplitude of matter fluctuations for all partially ionized scenarios will therefore lie in 
between their standard recombination (SR) and no recombination (NR) values as displayed in Fig. 2. All 
models are normalized to the COBE DMR detection at 10°. As is readily apparent from this figure, the value 
of 08, the present day mass overdensity on the 8h~ x Mpc scale, fits a function of the form cr&(n) = ae hn where 
the constants a and b are given for standard and no recombination in Tab. 1. The value of as for a specific 
n must fall within the range given by Fig. 2. Moreover, in an unbiased scenario of galaxy formation, as 1. 
Numerical simulations (Cen, Ostriker, & Peebles 1993) indeed seem to indicate that bias at 8/i _1 Mpc is low, 
b 1.1. Thus, for 0.1 < fio < 0.2 and 0.5 < h < 1.0, this implies that the spectral index must lie within the 
range —1.2 < n < —0.3 for any ionization history. The constraint for a specific ionization history is of course 
much more stringent. Growth during the period where the universe is transparent can be incorporated in a 
function which is nearly independent of the index n, 

a 8 (zi, k, a i ^ R) F(z i , Xi). (1) 



We will give the numerical values of this function for various fio and h in §3. A small n dependence 
arises at the SR extreme due to the change in the COBE DMR normalization with ionization history. 
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TAB. 1. Fitting parameters for as(n) = ae hn for standard recombination (SR) and no recombination (NR). 
These values can be used in conjunction with Fig. 7 to obtain as for partially ionized models. In 
compact object dominated models, as is nearly the SR value. 

2.2 Radiation Anisotropics 

Now let us consider the radiation power spectrum. We will separate the contributions to the final 
anisotropy into primary fluctuations, generated at standard recombination, and secondary fluctuations, gen- 
erated afterwards. Doppler shifts from scattering off electrons in small scale bulk motion will create ex- 
ceedingly large temperature inhomogenieties on the scattering surface at recombination, due to the large 
amount of small scale power in these models. As the radiation then free streams, the inhomogenieties become 
anisotropics under the horizon scale and give rise to primary temperature anisotropics on scales of several 
arcminute to degrees today. When the universe reionizes at z,, Thomson scattering off electrons isotropizes 
the photon distribution and reduces these small scale primary fluctuations to 
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is the optical depth due to Thomson scattering, where n e is the electron number density and a? is the 
Thomson cross section. For scales above the horizon, the fluctuations are still inhomogeneities, and being 
isotropic, do not damp due to rescattering. Requiring that the anisotropics satisfy the upper limit on 
arcminute scales from the OVRO experiment (AT/T)ovro < 2.1 x 10 -5 (Readhead et al. 1989), we set a 
lower bound on the optical depth T m , n for any given model. In Fig. 3, we plot this minimal optical depth. 
Notice that it in all cases T m , n is of 0(1) and is only weakly dependent on fio- 




FIG. 3. Primary constraint from OVRO as a function of spectral index n. Thick lines represent h = 1.0 
whereas thin lines denote h = 0.5 for two different values of flo> = 0-1 (dashed) and fio = 0.2 (solid). 
There is a minimal optical depth necessary to damp primary anisotropics from standard recombination 
to levels under the OVRO upper limit of (AT/T)ovro < 2.1 x 10 -5 . The minimal optical depth is 
only a weak function of flo- 

Up until now we have neglected secondary anisotropics which are generated after standard recombina- 
tion. Many effects contribute to the final anisotropy: the Doppler effect on the new last scattering surface, 
the Vishniac (second order Doppler) effect (Ostriker & Vishniac 1986, Vishniac 1987), the adiabatic growth 
of intrinsic photon fluctuations, the ordinary Sachs- Wolfe effect on the new last scattering surface, the en- 
tropy effect, and the integrated Sachs- Wolfe effect (both the curvature and entropy components) (Sachs & 
Wolfe 1967). The last four effects can be treated in a simple and unified way by examining the Boltzmann 
equation for the photons (Hu & Sugiyama 1994). We will refer to the combination of these four effects as 
the total Sachs- Wolfe effect. 

For the Doppler type effects, the amplitude and shape of the spectrum depend sensitively on when last 
scattering occured, which is itself determined by x,. On scales smaller than the horizon at last scattering, 
the Doppler effect is damped due to the cancellation of redshifts and blueshifts from photons that scattered 
from the front and back of a perturbation. Since the horizon grows with time, this damping effect, due to 
the finite thickness of the last scattering surface, is maximized by having the most recent last scattering, 
i.e., the maximal x,. In Fig. 4, we plot the effects of varying x, on the temperature anisotropics. Note that 
the CV's are related to the ensemble average of the temperature fluctuations as 

where Wi is the experimental window function taken from White, Scott, & Silk (1994). In general, increasing 
x, moves the peak of the temperature distortions to larger scales and decreases the amplitude. There is 
a small effect that can increase the amplitude however. In a highly ionized universe, the drag epoch is 
somewhat earlier than the last scattering epoch, defined as the time when optical depth reaches unity: 
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Zdrag > %u ~ 30(f2o/0. 1) 1 ^ 3 (0.05/a;j'f2/GM^) 2 ^ 3 • The growth of velocities between the two epochs can lead 
to slightly larger temperature fluctuations due to the Doppler effect. As one can see from Fig. 4, this is a 
very small effect. 




FIG. 4. Dependence of temperature anisotropics on the ionization fraction x, after the reionization epoch 
Zi = 500. Fluctuations are damped more severely and on larger scales in high x, models due to the 
increasing thickness of the last scattering surface. The peak of fluctuations correspondingly moves to 
larger scales. Aside from a small effect due to growth between the drag epoch and last scattering, 
raising x, will always minimize fluctuations. 

Since the first order effect is strongly damped at sub-horizon scales, second order effects can play a role 
(Ostriker & Vishniac 1986, Vishniac 1987). Mode coupling between the density perturbations and the bulk 
velocity, called the Vishniac effect, can give rise to a significant contribution on arcminute scales. Moreover, 
it has been shown that this is the dominant second order effect (Hu, Scott, & Silk 1994). However since it 
is proportional to the square of the density fluctuations, it is only important if perturbations have grown to 
a significant level by last scattering. This implies that the second order effect will decrease rapidly as one 
lowers x, to make the last scattering earlier. We shall quantify this constraint in §3 by using the techniques 
developed by Efstathiou (1988) and extended by Hu, Scott, & Silk (1994) and Chiba, Sugiyama, & Suto 
(1994) to open universes. 

The integrated part of the total Sachs- Wolfe effect (Sachs & Wolfe 1967), which occurs since the photons 
are traveling through a time dependent potential, can also play a role at small scales for fully ionized models 
(Hu & Sugiyama 1993). The relative importance of the small scale effect decreases in partially ionized 
universes due to an increase in the Doppler term. On the other hand, at large scales the total Sachs- Wolfe 
effect is the dominant term in all cases. Since such fluctuations are always above the Jeans length, this effect 
is independent of the ionization history. Thus, the low multipole moments of the temperature anisotropy are 
nearly identical for all models. This then implies that the COBE DMR normalization is essentially fixed. It 
is also interesting to note that the slope of the spectrum at large scales is only weakly dependent on fio and n 
(Sugiyama & Silk 1994). Compared with the spectral index of flat adiabatic models n a d, low fio isocurvature 
models predict n ad ~ 2 in contrast with the inflationary prediction of n a d ~ 1. Recent indications of a steep 
adiabatic slope in the COBE DMR maps of n ad = 1.59t^5 (Bennett et al. 1994) may argue in favor of this 
feature. 

On intermediate scales, i.e., scales near the horizon at last scattering, the adiabatic growth of the 
coupled photon-baryon system plays a role. As mentioned above, a curvature perturbation which grows as 
the adiabatic mode will dominate as the fluctuation comes within the horizon. Therefore, the monopole 
component of the photon distribution, i.e., the energy density fluctuation in the photons, will grow with the 
baryons if the universe is ionized. After last scattering, this perturbation is transferred to anisotropics at 
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I ~ 100. Hence at intermediate scales, the photon fluctuations will be the largest for early ionization. The 
opposite is true for the small scales at which the Doppler effect dominates. For early ionization, the baryon 
velocity cannot grow due to Compton drag. Therefore at smaller scales, fluctuations will be lower for early 
ionization. In Fig. 5 we display the dependence of the anisotropics on z,. We have fixed x, = 0.1 to isolate 
this effect. In all these models, the last scattering surface consequently has the same redshift and thickness. 
As one can see, intermediate scale fluctuations are minimized by having low z,. 




FIG. 5. Dependence of temperature anisotropics on the reionization epoch z, with fixed last scattering 
surface (x, = 0.1 z\ s k, 110). In models with high z,, the photon fluctuations have had longer to grow 
adiabatically with the matter on large scales (£ < 100). On small scales, velocities are prevented from 
growing by Compton drag, and thus the Doppler effect is smaller for high z, models. 



2.3 Spectral Distortions 

Ionization also produces spectral distortions by the SZ effect (Zel'dovich & Sunyaev 1969) in all rea- 
sonable scenarios, since the intergalactic medium must have been heated by some process such as collisional 
ionization. A conservative estimate of the necessary temperature is T e > 5000K (e.g. Cen, Ostriker, & 
Peebles 1993 find T e « 10000 - 20000K). As CMB photons scatter off this hot medium, a spectral distortion 
of the Compton-t/ type will be established as low frequency photons are shifted to higher frequencies. The 
Compton-t/ parameter is given by 

,= /i!^f* (4) 

where the overdot is a time derivative, T = To(l + z), and To = 2.726 ± 0.005 is the present temperature of 
the CMB (Mather et al. 1994). For constant T e , the integral is approximately, 

««8.4x 10~ 7 r( Te 1 (1 - - — Zi) 

V 50001*7 V 5T e 7 (5) 

< 2.5 x 10 -5 , 

where the constraint comes from the COBE FIRAS experiment (Mather et al. 1994). For T e = 5000K, the 
FIRAS limit becomes r < 30/(1 — 0.003zj). Combined with the minimal optical depth needed to suppress 
primary fluctuations (see equation (2)), the ionization history is severely constrained. 
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2.4 Compact Objects 

Now let us consider the effect of having a significant fraction of the universe in compact objects that 
behave essentially as cold dark matter. The clustering of compact objects is independent of the ionization 
history. Thus, if the compact objects dominate dynamically over baryons in the intergalactic medium, i.e. 
fic ^> &igm, the baryons would have fallen into the compact object wells after the drag epoch. If the 
universe is transparent between recombination and compact object formation, the transfer function will be 
nearly identical to the standard recombination one. Note that the power spectrum in these models will have 
prominent oscillations on small scales. The peak will also be at smaller scales than for the partially ionized 
models. Prior treatments have on the contrary assumed that the power spectrum is identical to the fully or 
partially ionized spectrum (Cen, Ostriker, & Peebles 1993). In these models, as also takes on the value of 
the standard recombination result displayed in Fig. 2. 



- fi =0.2 h = 0.5 




FIG. 6. Dependence of temperature anisotropics on compact objects fic = fio — fi/GM- In all models, 
compact objects form at z c = 800, and the universe is transparent until z, = 400. For a fixed optical 
depth Xi^liGM /&o = 0.045, raising fic makes drag effects on the baryons larger. On large scales the 
photon "velocity" (i.e. dipole) is larger than the baryon velocity. The converse is true on small scales. 
Thus the baryon velocity will be kicked higher at large scales and dragged lower on small scales. The 
anisotropics suffer a similar effect due to feedback through the Doppler effect. 

The behavior of the radiation fluctuations is more complicated however. First of all, lowering the density 
of matter in the intergalactic medium makes last scattering earlier thereby creating larger fluctuations at 
small scales (see Fig. 4). To isolate the additional effects of compact objects, we will fix the optical depth 
t(z) in these models by requiring that x^igm and z, be the same for all models (see Fig. 6). This fixes 
the location of the last scattering surface and the suppression of the primary anisotropics. However, the 
Compton drag on any given baryon is independent of fi/GM and is larger for models with high x,. If the 
baryon velocity is greater than the photon "velocity" (dipole), then Compton drag will decrease the velocity. 
On the other hand, if the photon velocity is larger, the drag will turn into a kick and increase it. 

Before recombination, the baryon and photon velocities are comparable due to tight coupling. During 
the free streaming epoch between recombination and reionization, on scales smaller than the horizon at 
recombination, the matter perturbations have already joined into the adiabatic growing mode and will just 
grow in linear theory. The photon dipole however damps due to free streaming. Thus, at small scales the 
photon velocity at reionization is smaller than the baryon velocity. On larger scales, the baryon velocity 
decays slightly after recombination as it attempts to join the growing mode of adiabatic perturbations. Thus, 
the photon velocity at large scales is larger than the baryon velocity. Since last scattering happens relatively 
near the drag epoch in partially ionized models, there is a corresponding feedback to the temperature 
anisotropics due to the Doppler effect. Models with lower fi/GM (higher fic) will therefore experience larger 
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fluctuations at large scales I < 100 due to the kick imparted by the photons and smaller fluctuations at 
small scales due to the drag associated with the photons. Notice that in Fig. 6, there is also a contribution 
on very small scales I > 1000 due to the primary anisotropy. All models show the same amount of primary 
fluctuations since we have fixed the optical depth (see equation (2)). To minimize fluctuations at I < 100, 
drag must be minimized i.e., for fixed optical depth we want no compact objects. 

In summary, raising x, will uniformly decrease anisotropics, whereas lowering z, will decrease anisotropics 
at the intermediate (£ ~ 100) scale and increase them at the small scale. Furthermore, x, and z, must be 
low enough to satisfy constraints on spectral distortions (see eqn. (4)) but high enough to damp the primary 
fluctuations (see eqn. (2)). Since these dual constraints limit the range of possible optical depths, the addition 
of compact objects will have the added effect of raising intermediate scale fluctuations and lowering small scale 
fluctuations. Furthermore, since minimal fluctuations are generated by the latest allowable last scattering, 
all compact object models have intermediate scale temperature anisotropics larger than this minimal model. 
The relative complexity of this picture is just a function of the number of free parameters in this theory and 
unfortunately cannot be avoided. 

3. CONSTRAINTS AND PREDICTIONS 




FIG. 7. Constraint plane for the ionization history with n = —0.5. The "primary" constraint arises from 
requiring the standard recombination anisotropics be damped sufficiently to escape the OVRO limit 
of (AT/T)ovro < 2.1 x 10 -5 on the scale of several arcminutes (see equation (2) and Fig. 3). The 
Compton-t/ constraint assumes a conservative minimal temperature for ionization of T = 5000K. The 



10 



Vishniac constraint from the ATCA experiment (AT/T)atca < 0.9 x 10 -5 on arcminute scales con- 
strains models with recent last scattering and high normalization. Dashed lines represent the constraints 
if the upper limits of the corresponding experiment are improved by a factor of 2. Approximate contours 
of constant as are given (thin solid lines, with as = 1 bolded), which also define the ionization function 
F(zi, x,) (see equation (1) and text). The value of as for any n can be generated with the aid of Fig. 2 
or Tab. 1. Minimal fluctuations at I < 100 are given by the model with largest allowable x, and lowest 
allowable z,. The models of Fig. 9, which are essentially these minimal models, are marked out with an 
asterisk on this figure. Corresponding temperature anisotropics are given in Fig. 9 and Tab. 2. 

We will now quantify the considerations of the previous section in order to constrain the ionization 
history, spectral index, and the formation of compact objects and make model predictions for various degree 
scale anisotropy experiments. As discussed above, the constraint from primary fluctuations and the Compton- 
y parameter work in opposite senses: the former requires high optical depth, the latter low. The Vishniac 
effect constrains models with a high amplitude of present fluctuation (large as) and late last scattering. The 
ATCA measurement on arcminute scales of (AT/T)atca < 0.9 x 10 -5 (Subrahmanyan et al. 1993) severely 
constrains high fio and h models through this second order effect (Hu, Scott, & Silk 1994). 




FIG. 8. Constraint plane for the ionization history with n = —1.0 (see Fig. 7). 

Let us first consider the case of no compact objects. In Fig. 7, we combine the above constraints for four 
models which bracket the interesting values of fio and h. These models all have the favored spectral index of 
n = —0.5 (Cen, Ostriker, & Peebles 1993). We have also plotted the values of as as a function of ionization 
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history. Note that the value of as for any spectral index can be read from this figure by using Fig. 7 and 
equation (1). For example for fio = 0.1 and h = 0.5, F(z, = 500, x, = 0.1) = <r&(n = — 0.5, z, = 500, x, = 
0.1)/(jg NR \n = —0.5) = 2.6. Therefore <r&(n = —1.0, z, = 500, x, = 0.1) ps 0.45; whereas the actual value for 
this model is as = 0.46. Notice that if we require as ps 1, the FIRAS Compton-t/ constraint excludes nearly 
all models in this range of fio and h. We have also plotted in dashed lines the corresponding constraints if 
the upper limits on the three relevant experiments, OVRO, ATCA, and FIRAS are brought down by a factor 
of 2. In the case of the FIRAS Compton-t/ constraint, this dashed curve is also the constraint for the more 
realistic value of T e = 10000K in equation (4). Viable models therefore must have a steeper spectrum than 
n ps —0.5. Interestingly enough, numerical simulations show that n ps — 1 fully ionized models are consistent 
with large scale structure formation (Suginohara & Suto 1992). In Fig. 8, we plot the same constraints for 
n = —1.0, and we see that reasonable values of as fall within the allowed region. 
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TAB. 2. Comparison of predicted anisotropics for partially ionized models, (AT/T) rms in units of 10 -5 , with 
observational detections. Bold face numbers show predictions which are in conflict with observation at 
the 95% CL. The first 8 entries represent models with "minimal" fluctuations at £ ~ 100 of Fig 9. The 
absolute minimal fluctuations can be slightly smaller (~ 10%) due to growth between the drag epoch 
and last scattering (see text). All models are therefore excluded by the SP91 observation at the 95% 
limit. Models with high h and n are also inconsistent with most other experiments. Latter entries are 
chosen to be consistent with as ~ 1- A standard CDM with h = 0.5 and = 0.06 is also shown for 
comparison. For the experimental values, TENE is the quoted rms value; ARGO, MSAM2, MSAM3 
are scaled gaussian autocorrelation function results; SP91, PYTH, MAX are effective rms values for 
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an analysis that properly takes into account correlations between data points for a Harrison-Zel'dovich 
spectrum (White, Scott, Silk 1994). 

Now let us consider the radiation fluctuations on degree scales. We have shown in section §2 that the 
way to minimize fluctuations on the I < 100 scale is to maximize x,, minimize z, and have no compact 
objects. In Fig. 9, we plot the minimal fluctuations for the eight models of Figs 6 and 7. These models are 
all'm conflict with the SP91 13pt scan of Schuster et al. (1993) at greater than the 95% CL (see Table 2) and 
have significantly larger intermediate scale anisotropics than the standard CDM model. Since these models 
have the minimal fluctuations at the SP91 scale of I ~ 70, we conclude that if the results of this experiment 
are confirmed, all models regardless of ionization or compact object formation are excluded. 




FIG. 9. Minimal temperature anisotropics at I < 100. All models are inconsistent with the upper limit from 
SP91 (see Tab. 2, Schuster et al. 1993). Therefore all thermal histories are also inconsistent, including 
those with compact objects. Fluctuations are also significantly larger than in the standard CDM model 
(thick dashed line, h = 0.5, = 0.06 ) even in these minimal models. 




FIG. 10. Temperature anisotropy for partially ionized models with <r% ps 1. Models with early last scat- 
tering x, < 0.01 overproduce fluctuations since the last scattering surface is nearly at the standard 
recombination epoch. Corresponding predictions for various experiments are given in Tab. 2. 
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Unfortunately, agreement between degree scale experiments is marginal at best. In particular, PYTH 
which is at a very similar angular scale as SP91, detects a somewhat larger fluctuation. Furthermore, MAX 
GUM (Gundersen et al. 1993) has detected significantly larger fluctuations on a slightly smaller scale. In 
Table 2, we also compare the predictions of these models with several other experiments. Only h = 1.0 
and n = —0.5 models are excluded at the 95% CL by the PYTH experiment and none are excluded by 
MAX GUM. In fact, the MAX GUM and TENE results favor isocurvature baryon models over standard 
CDM. However, only h = 0.5 and n = —1.0 models satisfy ARGO and MSAM2. Given the current state of 
confusion, it is perhaps too early to declare isocurvature baryon models excluded on the basis of the SP91 
13pt scan (or any individual experiment). 

We therefore also show the predictions for realistic models which fit the requirement that <r& 1 (see 
Fig. 10 and Tab. 2). Models with low ionization fraction x, < 0.01 tend to overproduce fluctuations on 
sub-degree scales. Although the optical depth is high enough to erase primary fluctuations, in these models 
the last scattering surface is so close to z = 1000 that equally large fluctuations are generated on the new 
last scattering surface. For models with a more recent last scattering surface, we obtain values closer to the 
minimal ones described above. 

1000 



800 
600 
400 



200 



0.01 0.1 1 

x Agm/^o 

FIG. 11. Example of the modification of the constraint plane (Figs. 7 & 8) to the compact object dominated 
case. The primary and Compton-t/ constraints come from optical depth and thus x, — > x,£Iigm /&o- 
The Vishniac effect depends additionally on the amplitude of the fluctuations, i.e. <r%. Since as is fixed 
to be near the standard recombination value in these models, the Vishniac constraint is the same for a 
fixed x,£Iigm The other panels of Figs. 7 & 8 can be generated in a similar way. 

In models which are dominated by compact objects, as is nearly independent of the ionization history 
since the baryons fall into the compact object wells by the present. The values for as are the same as 
the standard recombination ones of Fig. 2. Moreover, the constraint diagrams (Figs. 7 and 8) are also 
easily generalized to the compact object dominated case. Since the primary fluctuation suppression and the 
Compton-t/ parameter are purely functions of optical depth, one only needs to replace x, with x,£Iigm /&o- 
The Vishniac effect behaves similarly except that it depends sensitively on the magnitude of the fluctuations 
at last scattering. It is therefore determined from the value of as, which is fixed in these models, and Xi 
which determines the epoch of last scattering. An example of a constraint plane is displayed in Fig 11. The 
other planes may be generated in a similar fashion. 

By requiring consistency with as = 1, we can fix the spectral index n in compact object dominated 
universes. We will now take the more general scenario where after compact object formation at z c k, 800, the 
universe can be partially ionized, e.g., to x c = 0.1 (Gnedin & Ostriker 1992, Cen, Ostriker, & Peebles 1993). 
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At the galaxy formation epoch z,, the universe becomes fully ionized x, = 1.0. Again, the constraints can be 
easily generalized to this model by noting that the primary and Compton-t/ limits are functions of optical 
depth. We take z, = 5, 20, 50, 800 as representative cases. The last model corresponds to a universe which 
was fully ionized at the epoch of compact object formation. In Fig. 12, we plot the anisotropy spectrum; the 
corresponding predictions for various experiments can be found in Table 3. The models with low secondary 
ionization redshift z, overproduce sub-degree scale fluctuations because, either the optical depth is too low 
to satisfy the constraint on primary anisotropics, or the last scattering surface is so early as to be nearly 
identical with standard recombination. If x c < 0.1, only models with z, > 100 seem to be viable, in conflict 
with numerical work which shows z, ps 20 — 30 in these models (Cen, Ostriker, & Peebles 1993). Notice that 
these compact object dominated models are nearly always inconsistent with PYTH, ARGO, and MSAM as 
well as SP91. Should these detections be confirmed, only models with lower small scale fluctuations, i.e., 
smaller n, will survive. Since in order to satisfy as ps 1, n must decrease with flo> this implies that higher 
fio values would be desirable. It however remains to be seen whether these low n models are consistent with 
large scale structure. 
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TAB. 3. Comparison of predicted anisotropics for compact object dominated models, (AT/T) rms in units 
of 10 -5 , which have n chosen to be consistent with as ~ 1 (see also Fig. 12). Bold faced values show 
predictions which are in conflict with observations at the 95% CL. These models have compact object 
formation at z c = 800 followed by partial ionization x c = 0.1 with Sl/au/fio = 0.9. All models are 
inconsistent with SP91 at greater than 95% CL. Models with low ionization redshift z, overproduce 
fluctuations due to thinness of the last scattering surface and lack of optical depth. Most models are 
inconsistent with a majority of degree scale experiments. Higher fio models are more consistent with 
observations. 
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FIG. 12. Temperature anisotropy for compact object dominated models fic/^o = 0.9 with as 1. The 
ionization level after compact object formation at z c = 800 is x c = 0.1 until z, at which point the 
ionization rises to x, = 1.0. Models with low ionization epoch z, tend to overproduce fluctuations since 
they either do not satisfy the optical depth constraint from primary fluctuations r < T m , n or have last 
scattering so early as to regenerate fluctuations as large as standard recombination ones. Predictions 
for various experiments are given in Tab. 3. 

4. DISCUSSION 

The CMB potentially provides the most serious constraint on the isocurvature baryon model. As we have 
seen, the COBE DMR detection and galaxy clustering on the 8/i _1 Mpc scale fixes the primordial spectral 
index for any given ionization history. Furthermore, the ionization history is no longer a free function. It is 
strongly constrained by spectral distortions, arcminute scale primary anisotropics, and the Vishniac effect. 
Given these constraints, degree scale anisotropics are the Achilles heel of this model. These anisotropics 
are significantly larger than their CDM counterparts and are already in conflict with some degree scale 
experiments. However, indications of a steep slope by COBE and the high detections of TENE and MAX 
argue in favor of this model. Even if the latter indications are confirmed, it remains to be seen if a fully 
consistent model can be built. Observational evidence of a smooth, featureless power spectrum (Peacock & 
Dodds 1994) may conflict with this model. This is especially true for compact object dominated universes 
where the peak of the matter power spectrum and its small scale oscillatory structure are prominent. Such 
models can be made more consistent with observations by taking compact object formation to be relatively 
late z c pa 100 and assuming yet another population of early forming objects which reionize the universe 
between standard recombination and compact object formation. The matter power spectrum in these models 
will be similar to our partially ionized ones. Of course, all this fine tuning of the ionization history in order to 
satisfy microwave background and large scale structure constraints should also arise naturally from physical 
processes in the model. Nevertheless, relatively large intermediate scale anisotropics is a robust prediction 
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of the model. With the expected rapid improvement in these measurements, we shall soon be able to 
make a definitive statement about this model regardless of uncertainties in the ionization history and other 
parameters. 
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